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Contributions to the neutron's EDM, d n , are calculated using a well-constrained Ansatz for 
the nucleon's Poincare covariant Fadde'ev amplitude. The momentum-dependent quark dressing 
amplifies the contribution from the current-quarks' EDMs; and dressed-quark confinement and 
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using the current-quark EDMs generated by a minimal three Higgs doublet model of spontaneous 
CP violation is close to the current experimental upper bound. 
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I. INTRODUCTION 

The action for any local quantum field theory is invari- 
ant under the transformation generated by the antiuni- 
tary operator CPT, which is the product of the inver- 
sions: C - charge conjugation; P - parity transformation; 
and T - time reversal. The combined CPT transforma- 
tion provides a rigorous correspondence between parti- 
cles and antiparticles, and it relates the S- matrix for any 
given process to its inverse process, where all the spins are 
flipped and the particles replaced by their antiparticles. 
Lorentz and CPT symmetry together have many conse- 
quences, among them, that the mass and total width of 
any particle are identical to those of its antiparticle. 

The decay of the CP-odd eigenstate K\ into a CP-even 
2ir final state demonstrates that the product of only C 
and P is not a good symmetry of the standard model. 
This entails that time reversal invariance must also be vi- 
olated and that too has been observed in detailed studies 
of the neutral kaon system [3] . The separate violation of 
CP and T invariance can be accommodated in the six- 
quark-flavour standard model: only five of the six pos- 
sible phases in the 3x3 Cabibbo-Kobayashi-Maskawa 
(CKM) matrix can be eliminated via phase rotations of 
the quark flavours - rotating all the quarks through the 
same phase cannot affect the CKM matrix - leaving one 
CP violating phase. (In the four flavour theory, by con- 
trast, the 2x2 CKM matrix has three possible phases, all 
of which can be eliminated.) However, some aspects of 
the standard model are not completely satisfactory, no- 
tably the understanding of the magnitude of the direct 
CP violating parameter e' j|, ||, ||] . 

Identifying phenomena that are inexplicable in the 
standard model is an important focus of contemporary 
nuclear and particle physics, and while the neutral kaon 
system is the archetype for CP and T violation, the stan- 
dard CKM model predicts much larger effects for B- 



mesons || . A primary goal of the B- factories now being 
developed is to check these predictions. This is a new 
domain for testing the standard model. 

However, it has long been known that the possession 
of an electric dipole moment (EDM) by a spin i parti- 
cle would signal the violation of time-reversal invariance. 
(The existence of a dipole moment signals a spherically 
asymmetric distribution of charge.) Any such effect is 
likely small, given the observed magnitude of CP and T 
violation in the neutral kaon system, and this makes neu- 
tral particles the obvious subject for experiments: the ex- 
istence of an electric monopole charge would overwhelm 
most signals of the dipole strength. It is therefore natu- 
ral to focus on the neutron, which is the simplest spin-i 
neutral system in nature, and attempts to measure the 
neutron's EDM have a long history || JjJ. 

Very accurate measurements are made possible by the 
ability to produce and store ultracold 10 eV) polarised 
neutrons. Immersing such a sample in uniform magnetic 
and electric fields, the existence of a nonzero neutron 
EDM, d n ^ 0, manifests itself through a difference be- 
tween the precession frequency of the neutron's spin when 
the magnetic and electric fields are parallel as compared 
to when they are antiparallel. Contemporary applica- 
tions of this technique yield an upper bound 

\dn\ < 6.3 x l(T 26 ecm (90%C.L.). (1) 

(NB. e/(2M„) = 1.0 x 10" 14 ecm. Therefore, writing 
d n = eh n /(2M), where M is the neutron's mass and h n 
is its "gyroelectric ratio," Eq. ([[]) corresponds to \h n \ < 
6.0 xlO" 12 .) 

Equation (|l|) has proven to be an effective constraint. 
For example, the most general Lagrangian density for a 
local, renormalisable colour-S'C/(3) gauge theory would 
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contain a term: 



pa pa 
fjbi/ per ' 



(2) 



where F^ v is the non-Abelian field strength tensor and 9 
is an undetermined constant - the ll 9 angle." However, 
this term violates CP invariance and hence generates a 
nonzero neutron EDM. Its strength; i.e., the value of 6, 
is therefore constrained by Eq. (|l|), which yields a very 
low upper bound || ||, [l(J 



\e\ < (i - io) x 10" 



(3) 



Currently there is no satisfactory explanation of why this 
term is so strongly suppressed and that is the basis of the 
so-called "strong CP problem:" the goal is to find a rea- 
son why 6 should naturally be identically zero. (In the 
absence of topologically nontrivial gauge field configura- 
tions, Eq. (^) cannot contribute to the action: it is a 
surface term.) 

Herein we shall assume 9 = so that the phase in the 
CKM matrix is the only source of CP and T violation 
in the standard model. In this case a nonzero neutron 
EDM is proportional to the CKM phase. (Considerations 
relevant for 8 ^ are explored, e.g., in Refs. |3, [k|.) 

Over many years the experimental upper bound on 
d n has steadily decreased and its small value has also 
proven very effective in ruling out candidates for theories 
that enlarge the standard model. To illustrate how, we 
suppose for the moment that the neutron is a collection 
of three valence-quarks described by a symmetric SU (6) 
spin-flavour wave function. Then, by analogy with the 
magnetic moment, 



(4) 



where d u ,d are valence-quark EDMs. In the standard 
model the first nonzero contribution to a free quark's 
EDM appears at third order and involves a gluon ra- 
diative correction (i.e., 0(a s G 2 F ), for the same reason 
that flavour-changing neutral currents are suppressed: 
the GIM mechanism) so that §, |l| 



d u <d d < 10- J4 ecm 



(5) 



Using this in Eq. (Q) gives a result seven orders-of- 
magnitude less than the current experimental bound. 
This result is characteristic: other plausible mechanisms 
within the standard model, such as hadronic loop correc- 
tions, also give a very small value. However, the standard 
model is peculiar in this regard and candidates for its ex- 
tension typically contain many more possibilities for CP 
and T violation, which a priori are not constrained to be 
small. Thus Eq. (Q) is an important and direct constraint 
on these extensions because Eqs. (Q), (||) indicate that the 
standard model contribution to d n cannot possibly inter- 
fere at a level that could currently cause confusion. For 
example, as our calculation will show, the viability of the 



minimal model of spontaneous CP violation 12 , which 
involves three Higgs doublets, is endangered by Eq. (0). 

Our interest in the neutron's EDM stems from a desire 
to explore the validity of Eq. (||) in the sense that, irre- 
spective of the origin of the valence-quark EDMs, how 
are they related to the EDM of the bound state? This 
question has recently been explored H [l(], using QCD 
Sum Rules flli]] . Our analysis, however, will employ in- 
stead a recently developed, well-constrained Poincare co- 
variant, bound-state picture of the neutron jl^, 1C I [Tt], 



|l8| , [[9], |2fJ] and therefore affords a necessary complement. 

The nonzero neutron charge radius is a clear indica- 
tion that a symmetric SU(6) wave function is inadequate 
for the neutron and there are significant additional weak- 
nesses. For example, in making the connection between 
Eqs. (ji|) and (||) no consideration is given to the neces- 
sary momentum-dependence of the dressed-quark mass 
function, which is a longstanding prediction of Dyson- 
Schwinger equation (DSE) studies; e.g., Ref. (|l], p2| , 
that has recently been confirmed in lattice-QCD simula- 
tions |23[| , nor for confinement and the concomitant fea- 
ture that dressed-quarks in the neutron are not on shell. 
(Reference E1J provides an heuristic guide to DSEs and 
References p4L [2^] give an overview of their contempo- 
rary application, with particular emphasis on continuum 
strong QCD.) 

Hitherto, lacking a Poincare covariant bound state 
picture of the neutron, these aspects of dressed-quark 
behaviour and hadron compositeness have only been 
explored in studies of the dipole moments of the p- 
meson [^6[ |27| , a bound state for which sound DSE mod- 
els exist. Assuming an explicit EDM contribution to the 
photon-currcnt-quark vertex: 



2 m f 



(6) 



where: e is the positron charge; ?n/ =u ,d are the current- 
quark masses; e Ui d, their electric charge fractions; h Ut d, 
their gyroelectric ratios; and q^, the momentum transfer, 
an analogue of Eq. (||) for a pointlike p-meson is |26j 



d u + dj . 



(7) 



In this case the direct calculations showed that Eq. (Q) 
underestimates the EDM of a bound state p-meson, 
which is composed of a dressed-quark and a dressed- 
antiquark, and described by a Bethe-Salpeter amplitude, 
by as much as two orders-of-magnitude. That can mate- 
rially affect the inferred constraints on extensions of the 
standard model. 

A quantum field theoretical description of the nuc- 
leon as a composite of a quark and nonpointlike 
colour-antitriplet diquark was proposed and explored in 
Refs. [[[5], jig] ) and recent studies |l7]] have shown that 
this approach, based on a Poincare covariant Fadde'ev 
equation, is capable of providing a good description of 
the spectrum of octet and decuplet baryons. With this 
foundation, a product Ansatz for the nucleon's Fadde'ev 
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amplitude has been employed [L9| |2C]] to describe a 
wide range of elastic nuclcon form factors. This is the 
model we use to explore the implications of Refs. (^(| Q 
for the neutron's EDM. 

In Sec. H we describe our model in detail, including a 
discussion of the form of a CP and T violating coupling 
of a photon to a dressed-quark. An analogous dressed- 
quark-gluon coupling is also admissible Q and may yield 
an equally important contribution to d n We ne- 

glect it and hence ours is not a complete calculation of 
d n . Nevertheless, this and like terms are additive, and 
their omission does not qualitatively affect our discus- 
sion nor the points we wish to emphasise. Our results 
are presented in Sec. [II, and we report more than just 
the EDM so as to establish a context for that result and 
demonstrate the level of accuracy to be anticipated using 
our model. Section IV is an epilogue. 



II. MODEL ELEMENTS 



A. Dressed-quarks 



The general form of the dressed-quark propagator 



is ra 



S(p) 



-17 • pa v (p 2 ) + <J s (p 2 ) , 



(8) 
(9) 



and it is a model-independent result of quark-DSE stud- 
ies that the wave function renormalisation and dressed- 
quark mass: 



Z(p 2 ) = 1/A(p 2 ) , M(p 2 ) = B(p 2 )/A{p 2 ) , 



(10) 



respectively, exhibit significant momentum dependence 
for p 2 < 1 GeV 2 , which is nonperturbative in origin. (See, 
e.g., Ref. (^9|.) This behaviour is a longtime prediction 
of DSE studies j2l|] and has recently been observed in 
lattice-QCD simulations p3| . The infrared enhancement 
of M(p 2 ) is an essential consequence of dynamical chi- 
ral symmetry breaking (DCSB) and is the origin of the 
constituent-quark mass. With increasing p 2 the mass 
function evolves to reproduce the asymptotic behaviour 
familiar from perturbative analyses, and that behaviour 



is unambiguously evident for p 2 > 10 GeV 2 22 



While numerical solutions of the quark DSE are read- 
ily obtained, the utility of an algebraic form for S(p) 
when calculations require the evaluation of numerous 
multidimensional integrals is self-evident. With this 
in mind an efficacious parametrisation of the dressed- 
quark propagator, which exhibits the essential features 
described above, was introduced in Ref. Q and has 
been used extensively in studies of meson properties; e.g., 
Refs. iol P, M |3lL P, PL M. We use it herein. 



The parametrisation is expressed via 

a s (x) = 2mf(2(i + m 2 )) (11) 
+T(b 1 x) F(b 3 x) [b + b 2 T{ex)\ , 

a v (x) = [l-T{2{x + m 2 ))} , (12) 

x + m z 

with x = p 2 /X 2 1 to = to/A, J-(x) — [1 — exp(— x)]/x, 
as(x) = \as(p 2 ) and <Jv(x) = \ 2 <Jv{p 2 )- The mass- 
scale, A = 0.566 GeV, and parameter values 



to b bi b 2 



0.00897 0.131 2.90 0.603 0.185 



(13) 



were fixed in a least-squares fit to light-meson observ- 
ables plf. The dimensionless u — d current-quark mass 
in Eq. (|13|) corresponds to 



5.1 MeV. 



(14) 



(e = 10~ 4 in Eq. dTT| ) acts only to decouple the large- 
and intermediate-p^iomains.) 

Dynamical chiral symmetry breaking is expressed in 
the parametrisation. It gives a Euclidean constituent- 
quark mass 



0.33 GeV, 



(15) 



defined j?2| as the solution of p 2 — M 2 (p 2 ), whose 
magnitude is typical of that employed in constituent- 
quark models |37| and for which the value of the ratio: 
M^ d /m — 65, is characteristic and definitive of light- 
quarks |Q. In addition, DCSB is manifest through a 
nonzero vacuum quark condensate 



lGeV 2 



A' 5 



3 &o 
47r 2 bi 6j 



In 



lGeV 2 
A 1 

QCD 



(0.221 GeV) 3 , 
(16) 



where we have used Aqce> = 0.2 GeV. The condensate is 
calculated directly from its gauge invariant definition |3q | 
after making allowance for the fact that Eqs. (p"T|), ( |12] ) 
yield a chiral-limit quark mass function with anomalous 
dimension 7 m = 1. This omission of the additional 
ln(p 2 /Aq CD )-suppression that is characteristic of QCD 
is a practical but not necessary simplification. 

Motivated by model DSE studies |§, Eqs. (Jll]) , (|l|) 
express the dressed-quark propagator as an entire func- 
tion. Hence S(p) does not have a Lehmann representa- 
tion, which is a sufficient condition for confinement [^0[ . 
Employing an entire function for S(p), whose form is 
only constrained by the calculation of spacelike observ- 
ables, can lead to model artefacts when it is employed 
directly to calculate observables involving large timclikc 
momenta An improved parametrisation is there- 

fore being sought. Nevertheless, on the subdomain of 
the complex plane explored in the present calculation 
the integral support provided by an equally efficacious 
alternative cannot differ significantly from that of this 
parametrisation, which explains why we use it herein. 
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B. Nucleon 

The idea that diquark correlations play a significant 
role in baryon structure and interactions is almost as old 
as that of quarks themselves . It was a motivation for 
the meson-diquark bosonisation enunciated in Ref. J43| , 
which provides a picture of baryons as dressed-quark- 
diquark composites, and the subsequent derivation |l5| of 
an homogeneous, Poincare covariant Fadde'ev equation 
for baryons that exploits the role of diquark correlations. 

Our picture of the nucleon is based on the latter ap- 
proach. We represent the nucleon as a relativistic three- 
quark bound state, involving a nonpointlike, Lorentz- 
scalar diquark correlation, via a product Ansatz for the 
Fadde'ev amplitude: 



*s(pi; o>i, n) = e ClC2C3 A 0+ (K) 

x \r 0+ £ P[12] ;K)]™ 2 m-p)u{P)]Z 



(17) 



where: (ij ■ P + M)u(P) = 0, with P = pi + P2 + P3 —■ 
P{i23} the nucleon's total momentum and M its mass; 
£cic 2 c 3 is the Levi-Civita symbol that provides the colour- 
singlet factor; K = p x + p 2 =: P{i2}, P[i2] ■= Pi ~ P2, 
£ = (p{i2} — 2j>3)/3; and (aj, t^) are the quark spinor and 
isospin labels. 

Equation ( p"7| ) describes the general form of the ampli- 
tude in the scalar diquark subspace. In this equation 
A 0+ (K) is the pseudoparticle propagator for a scalar 
diquark formed from quarks 1 and 2, and T + is a 
Bethe-Salpeter-like amplitude describing their relative 
momentum correlation. As explained, e.g., in Ref. p4| , 
these functions can be obtained from an analysis of 
the quark-quark scattering matrix. However, follow- 
ing Refs. |ll| |l9|, p0| , for simplicity herein we employ 
parametrisations: 



A a+ (K 2 ) 



— T(K 2 /u; 2 Q+ ), 



(18) 



T Q+ (k;K) = Ci l5 ir 2 T(k 2 /u; 2 0+ ), (19) 
yv + 

where C = 7274 is the charge conjugation matrix and 
A/o+ is a calculated, canonical normalisation constant 
that ensures, e.g., that a (wd)-diquark has electric charge 



fraction (1/3) for K 2 



l o+- 



The parameters are a 



width, uj + , and a pseudoparticle mass, m + , which have 
ready physical interpretations: the length r + = 1/lu + 
is a measure of the mean separation between the quarks 
in the scalar diquark; and the distance l qq — l/m qq rep- 
resents the range over which a true diquark correlation 
in this channel can persist inside a baryon. (NB. The ab- 
sence of a particle-like singularity in the pseudoparticle 
propagator presented in Eq. ([l8|) is sufficient to ensure 
that the diquark is confined inside the baryon p0|.) 

The remaining element in Eq. (|l7|), -0, is a Bethe- 
Salpeter-like amplitude that describes the relative mo- 
mentum correlation between the dormant quark and the 
diquark's centre-of-momentum. Using Eqs. (jl^), (|l9|), ip 



can be obtained by solvinga Poincare covariant Fadde'ev 
equation for the nucleon p.q|: 



V>(fc; P) u(P) = -2 



d 4 e 



A (K e )T 0+ (k + £/2;K) 



(2tt)4 

xS(£ cx ) T T Q+ (£ + k/2;-K k ) 

xs(ii)ip(e-,p)u(p), 



(20) 



where K e = -£+(2/3)P, £ cx = -£-k-P/3, t x = £+P/3. 
For a positive energy nucleon, the solution has the general 
form 

m P) = h(£; P)l-±(i<y.l-l.p l) f 2 (£; P) , 

(21) 

where P 2 — — 1 and, in the nucleon rest frame, /i^ de- 
scribe, respectively, the upper/lower component of the 
dressed-nucleon spinor. 

To learn about tp we solved Eq. ( |20| ) and also its ex- 
tension to include an axial-vector diquark correlation, a 
sound foundation for the latter step having been laid in 
Ref. [0. We used the dressed-quark propagator de- 
scribed in Sec. p A| and, to further simplify and so ex- 
pedite the calculations, we retained only the first Cheby- 
shev moment of the functions /1.2 in Eq. (plf); i.e., we as- 
sumed fi, 2 (£; P) = hM e2 i P 2 )- In thi s way we found |(J 
that a simultaneous description of the nucleon and A 
masses is possible when the axial- vector diquark correla- 
tions are included. That description requires (in GeV) 



(22) 



where W]+ , are obvious analogues in the axial- vector 
channel of the scalar diquark parameters in Eqs. (|l8]), 
(|j~9|), and corresponds to (in fm) r + = 0.47, r x + = 0.23, 
/ + = 0.31, li+ = 0.18. (The last two columns in 
Eq. ( p2| ) are the calculated masses. A description of 
the A, for which J = 3/2, is obviously impossible with 
only a scalar diquark.) In our calculation the value of 
m + was taken from Refs. jl8|, [ljj and that of the ra- 
tio rnQ+fm^ = 0.78 from the Bethe-Salpeter equation 
studies of Ref. [^|, which is consistent with that ob- 
tained in lattice-QCD simulations pq |. The parameters 
w o+,i+ were then varied to fit -M/v,a- The observations: 
m +/m 1+ = 0.78 w 0.76 = M N /M A ; M^d + m o+ = °- 97 ; 
and M^ d +m 1+ = 1.19, are suggestive. However, our own 
and other [jl7[ ^7) Fadde'ev equation studies show that 
the nucleon contains a significant axial-vector diquark 
component (neglecting it, the nucleon mass is ~ 40% too 
large j20|) and hence the origin of the N and A masses 
and mass splitting does not lie simply in summing over 
constituent masses. 

For the purpose of developing an intuitive understand- 
ing, our eigenvector solution for the nucleon can ade- 
quately be approximated as [M 



w + 


m 0+ 


U>1 + 


m 1+ 


M N 


M A 


0.42 


0.64 


1.09 


0.86 


0.94 


1.23 



i-y-t-l-Pl) 



(23) 
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with calculated values of lo^ « 0.4 GeV (r^ = 0.49 fm) 



and R w 0.5. This makes clear that 



> 



r + 



> 



ri+, 



which is a necessary condition for the internal consis- 
tency of the quark+diquark Fadde'ev equation descrip- 
tion since it signifies that the mean separation between 
the quarks in the constituent diquarks is no more than 
the size of the nucleon. Furthermore, the value of R indi- 
cates that the lower component is a significant piece of a 
relativistic nucleon's spinor. (A/^ is a calculated normal- 
isation constant that ensures, e.g., that the proton has 
unit charge.) 

Following this study we can complete the specification 
of a well-informed product Ansatz. Equation (fT^), with 
Eqs. (|l8|), (|f9|), (p3|), provides a two-parameter model: we 
fix the values of R = 0.5 and m + w 0.64 GeV, motivated 
by the Fadde'ev equation studies, and allow lo + and ui^ 
to vary so as to obtain a least-squares fit to the proton's 
electric form factor, as described in Refs. |T^, [f9|, p0[ . 
With this expedient we have an algebraic scalar diquark 
model that provides accurate estimates of known observ- 
ables, as we show in Sec. HI, and hence can easily be used 



to obtain realistic constraints on the neutron's EDM that 
reflect the influence of those aspects of strong QCD that 
we identified in Sec. |j: DCSB; quark confinement; and 
hadron compositeness. 



C. Quark-Photon Coupling 

A calculation of the electromagnetic interaction of a 
composite particle cannot proceed without an under- 
standing of the coupling between the photon and the 
bound state's constituents. This is illustrated with par- 
ticular emphasis in Refs. |50| |33|, |4S[ ], which consider ef- 
fects associated with the Abclian anomaly. When quarks 
are dressed as described in Sec. II A, only a dressed- 
quark-photon vertex, T^, can satisfy the vector Ward- 
Takahashi identity: 



(24) 



where q = £\ — £ 2 is the photon momentum flowing 
into the vertex. The constraints that this identity and 
other features of a renormalisable quantum field theory 
place on the form of have been explored extensively 
in Refs. g§]. 

is the solution of an inhomogeneous BSE, and the 
pointwise behaviour of the solution has been elucidated in 
the numerical studies of Refs. [^o| . However, for our pur- 
poses we again prefer an efficacious algebraic parametri- 
sation and that is provided by |plf 



ii 7 -(4 + ^)A^,£!) + A s (^ 



A F (£l,e 2 2 ) 



F{l\)-F{li) 



(I 



(25) 
(26) 

(27) 



where F = A, B\ i.e., the scalar functions in Eq. (||). A 
feature of Eq. (|2^) is that is completely determined by 
the dressed-quark propagator. Furthermore, when using 
the dressed-quark propagator parametrisation, improve- 
ments on this vertex Ansatz can only modify our results 
by < 10 %, as illustrated, e.g., in Refs. [36, |5^| . 

Equations and Refs. @ consider only 

the CP preserving part of the dressed-vertex. When the 
possibility of CP and T violation is admitted, additional 
contributions are possible and the form that has most 
often been considered is that of Eq. (pi): 



p 2m 



(28) 



where: e is the positron charge, m is the u = d current- 
quark mass; Q — diag[2/3, — I /3] is the quark charge 
fraction isospin matrix; and H~ = di&g[h~(q 2 ), h^(q 2 )] 
defines an analogous quark EDM matrix, wherein we ac- 
knowledge that the quarks' EDMs may be momentum- 
dependent. Plainly, adding this term preserves Eq. (|2~i|). 

As discussed in Ref. J2?J, another contribution is pos- 
sible: 

r+(*i,*a) = 7^QH + (q 2 )l l5 
p 2 m 

x [(h + - & + 40 • q%\, (29) 

where: q 2 = 1; H + {q 2 ) is an obvious analogue of H~(q 2 )\ 
and adding this term also preserves Eq. (|24|). Note 
though that in any calculation where the quarks are as- 
sumed to be on-shell, and hence described by spinors for 
which (17 • I + m)u{£) — 0, then 

iu(£t) 75 (it + h)^ u{£ 2 ) = u(h) 75 a„u q v u{£ 2 ) , (30) 

and consequently the two structures are equivalent. How- 
ever, Eq. (30) is not satisfied by the dressed-quarks in the 
bound state nucleon: they are confined and may not even 
have a mass shell, so that in the general case the two con- 
tributions are distinguishable and must be treated sepa- 
rately. As an example, for the p-meson the different op- 
erator structures generate individual contributions that 
differ by < 20% from their average value p7j . 

In our calculations then we employ the following alge- 
braic Ansatz for the dressed-quark-photon vertex 



(£1 , £2 ) = Q r p (£1 , £ 2 ) + r+ (£, , £ 2 ) + r - {h , h ) 



(31) 



III. RESULTS 



The electromagnetic nucleon current is 
J^P',P) = ieu(P')A^q,P)u(P), (32) 
= ieu(P') ^F 1 {q 2 ) + ^- i a^q v F 2 {q 2 ) 

+ ^l^q,n(q 2 ))u(P), (33) 
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TABLE I: Calculated values of a range of well-known phys- 
ical observables. The "Obs." column reports experimental 
values [153 1 or values employed in a typical meson exchange 
model [154 1 . The remaining columns report our results, ob- 
tained using the Fadde'ev Ansatz parameters in Eq. (Bq) . 





Obs. 


Set 1 


Set 2 


(r v Y (fin*) 
(rnf (fm 2 ) 


(0.87) 2 
-(0.34) 2 


(0.78) 2 
-(0.40) 2 


(0.83) 2 
-(0.33) 2 


[UN) 


2.79 


2.83 


2.82 


[In (/LW) 


-1.91 


-1.61 


-1.62 


/^n / /^p 


-0.68 


-0.57 


-0.57 


gwNN 


13.4 


13.9 


15.4 


(tInn) (fm 2 ) 


(0.93-1.06) 2 


(0.63) 2 


(0.69) 2 




1.26 


0.98 


1.27 


<^> (fm 2 ) 


(0.68±0.12) 2 


(0.83) 2 


(0.77) 2 


QpNN 


6.4 


6.66 


8.60 


fpNN 


13.0 


15.8 


16.7 


guNN 


7-10.5 


12.2 


15.5 


fuiNN 




6.5 


5.2 



where: the spinors satisfy 7 ■ Pu(P) — iMu(P), u(P) 7 • 
P = iMu(P), with M = 0.94 GeV; R = P' + P and 
q ■ P = 0; and A M is the nucleon-photon vertex. F\ and 
F2 are the usual Dirac and Pauli electromagnetic form 
factors of the nucleon, in terms of which the electric and 
magnetic form factors are 



G E (q 2 



F^q 2 



4M 2 

G M (q 2 ) = F^q 2 ) + F 2 (q 2 ) 



(34) 
(35) 



The remaining term yields Ti. and hence the nucleoli's 
EDM form factor, which vanishes in the absence of a CP 
and T violating quark-photon coupling. Equation ( |33| ) is 
the general form because Eq. ([30j) is valid for an on-shcll 
nucleon. 

With the specification of the elements in Sec. ||, an im- 
pulse approximation calculation of the electromagnetic 
properties of the nucleon is straightforward using the ex- 
pression for A M given in the appendix. Each calculation 
requires the evaluation of a number of multidimensional 
integrals, which we accomplish using Monte-Carlo meth- 
ods, requiring a statistical accuracy of < 1 %. 

Our results for an illustrative range of well-known 
quantities are presented in Table |, which serves to ex- 
emplify the accuracy of our scalar diquark model of the 
nucleon. To display the sturdiness of the model we used 
two parameter sets (dimensioned quantities in GeV) 





R 


TO + 


w + 


lO-tp 


Set 1 


0.5 


0.62 


0.79 


0.23 


Set 2 


1.0 


0.63 


0.59 


0.18 



(36) 



with, in each case, the values of w + and tu^ determined 
via a least-squares fit to the proton's electric form factor. 
Set 1 is "preferred" because the value of R is fixed to 
that o btain ed in the Fadde'ev equation study described 
in Sec. [II B|. (The minor adjustment of m + is necessary 



to preserve neutrality of the neutron subsequent to the 
variation of u> +, aty, see next paragraph.) In compari- 
son with the form factors depicted in Refs. (l8[ |l^], both 
of which used R = and hence overlooked an impor- 
tant qualitative outcome of Fadde'ev equation studies, 
our calculated G r ^(q^) is much improved cf. the data, as 
seen, e.g., in Ref. EQ|. The other form factors are lit- 
tle affected. (NB. Only the electromagnetic properties 
reported in Table | can be calculated directly from the 
information provided herein. The remaining quantities 
were calculated as described in Ref. [|l9|: ttNN - using 
Eq (22) of that reference; vector meson - using Eq. (38); 
and axial form factor - using Eq. (60).) 

We remark that there is a quantitative discrepancy be- 
tween the values of uj + , ui^ obtained here and those ob- 
tained in the Fadde'ev equation study, Eqs. (p2|), (|2^). 
That is to be anticipated because here we are requiring 
a fit to wide range observables in a circumscribed model 
space; i.e., the discrepancy is an artefact of our scalar- 
diquark product Ansatz. However, our model's simplicity 
and illustrative utility outweigh the defect because in this 
and its earlier applications the qualitative implications of 
the defect are readily identifiable. 



A. Electric dipole moment 

Having established a context, we now report our results 
for the EDM. Using Eq. © one finds 



W = 5 dia g| E E «"?(») I' 



*=+.- 

\f=u,d 



"= + ,- 
f=u,d 



(37) 



an isospin matrix, where the calculated values of wj (N) 
express the dependence of the nucleon's EDM on its in- 
ternal structure. Clearly, because of isospin symmetry, 



«£(p)=<(n), <(p)=<(n) 



(38) 



The M/m multiplicative factor in Eq. ( p7| ) makes ex- 
plicit the importance of the DCSB mechanism in amplify- 
ing the contribution from dressed-quarks to the nucleoli's 
EDM. 

Existing experiments constrain Ti.(q 2 — 0) and in Ta- 
ble [n] we report the relevant values of Wj(n). (There 
are no cancell atio n s be tween contributions from differ- 
ent Alj Eqs. (A1)-(A7). The values are comparable in 



magnitude to their analogues obtained using QCD Sum 
Rules, e.g., Ref. fl3[|. ) The table makes clear the ex- 
tent to which Eq. (|30|) (the on-shell assumption) is vio- 
lated by the dressed-quarks in the nucleon bound state: 
for the doubly represented quark flavour the difference is 
quantitatively similar to that in the p- meson ^7|. The 
difference is much larger for the odd flavour quark. How- 
ever, that is likely an artefact of o nly r etaining a scalar 
diquark, in which case A^ in Eq. (A6) provides a sole, 
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TABLE II: Calculated values of the coefficients in Eq. (|^). 
(All quantities are dimensionless.) In the absence of dressed- 
quark confinement and off-shell effects, the entries in the two 
rightmost columns would be zero. The parameter values for 
Sets 1 and 2 are given in Eq. (pS6|). 

wj{n) w+{n) w d (n) w u (n) J_ ± " 1 



Set 1 0.560 
Set 2 0.647 



0.055 
0.043 



0.756 
0.878 



-0.111 
-0.102 



0.15 
0.15 



2.9 
2.5 



unmatched contribution to h . On the information avail- 
able we therefore judge that the magnitude of this effect 
is best estimated from the difference for the doubly rep- 
resented flavour. 

Using Eqs. @, @ and the Set 1 results from Table [| 
we obtain 



185.2- (-0.560 /i+ + 0.110 /i+ 



-0.756 hq - 0.222 ft-) , (39) 



hp = 185.2 - (-0.055/1+ + 0.111 h d 



1.12 h+ + 1.51 K 



(40) 



The ratio of the coefficients of hq and h~ in Eq. ( 139 ) is 
3.4, while the value of this ratio obtained from Eq. (Q) is 
2.0. The difference between these values gauges the ex- 
tent to which SU (6) spin-flavour symmetry is broken by 
quark + quark = diquark clustering in the nucleon bound 
state. This difference may change a little but will not 
vanish upon the inclusion of axial- vector diquark correla- 
tions in the nucleon's Fadde'ev amplitude. (The limited 
effect that axial-vector diquarks have on other observ- 
ables that do not involve cancellations between the A^- 
contributions |2(], |47[] makes us confident of this.) The 
h + contributions do not have an analogue perturbatively 
nor in the constituent quark model. 

We observe that hadronic loop insertions are unam- 
biguously an additive correction to our impulse ap- 
proximation and, based on analyses of other observ- 
ables ||^, ||^, ||^, ||^, |5?|, we expect that for realistic 
hadron masses their contribution will modify Eqs (39), 
(fh^) by < 10%. This suppression is due to the compos- 
iteness of the hadrons: away from the position of well- 
known kinematic singularities, the need to extrapolate 
bound state properties off-shell Q, for the evaluation 
of integrals over spacelike momenta; and the integration 
cutoffs applied by the hadrons' finite size, both lead to a 
quenching of the loop contributions. 

If we assume that CP-violating current-quark-level in- 
teractions yield hd 3> hy. as is the case, e.g., in Higgs 
boson exchange models [E7| , then 



h n = -40.6 (0.85 h+ + 1.15 hq) 



(41) 



and Eq. (|l|) applies the following bound on the e?-current- 
quark gyroelectric ratios: 



Using Eq. ( filf ) it is clear that, for hd 3> h u , d p < 0.1 d n . 

The DCSB mechanism, which is responsible for turning 
the current-quark mass i nto the constituent-quark mass, 
as reviewed in Sec. II A , is responsible for the suppres- 
sion of hd with respect to h n by the factor of ~ 40 in 
Eq. ([id]). This value is a lower bound on the magnitude 
of the effect. It will increase when using any reasonable 
Ansatz that incorporates strong interaction dressing in 
r+ analogous to that described by £p, in the CP 
and T preserving part of the vertex. We can estimate the 
scale of this effect by comparing u)~(p)(= u^(n)) and 
w d (p)(— u u ( n )) W1 th lattice estimates [|57| of the pro- 
ton's tensor charges: 6u = 0.839(60), Sd = -0.180(10). 
It is thus apparent that Eq. (^) overestimates the upper 
bound by less than a factor of two. 

The result in Eq. (JsJ> corresponds to \h^ + h^\ < 10~ 22 , 
so Eq. (^) does not challenge the standard model. How- 
ever, with j§ Re(e'/e) = (2.1 ± 0.46) x 10~ 3 , the calcu- 
lation in Sec. VI of Ref. p8| corresponds to a Weinberg 



model 1 12 prediction of [59 



|/i+ + /i d -| = (0.3~9.0)xl0- 



(43) 



Hence this model is threatened by our estimate of the 
bound, Eq. ([l2|), which incorporates a well-constrained 
modelling of the effects of DCSB, and the binding and 
confinement of dressed-quarks in the neutron. 

We emphasise that the results in Table [n| are indepen- 
dent of the means used to calculate hf. Hence Eq. ( J3S| ) 
can be applied to constrain any extension of the stan- 
dard model. Here we have only exemplified that using a 
particular candidate. 



IV. EPILOGUE 

The results are clear. Quantitatively - Our value for 
d n will be modified by the well-constrained inclusion of 
axial-vector diquark correlations. However, based on 
the effect this has on other observables, the correction 
should not exceed 15%. Qualitatively - 1) The scale 
of dynamical chiral symmetry breaking in QCD gener- 
ates a significant amplification of the contribution from 
a current-quark's EDM to that of the bound state con- 
taining it: the minimal enhancement factor is roughly the 
ratio of the constituent-quark and current-quark masses, 
and is understandable and calculable through the neces- 
sary momentum-dependence of the dressed-quark mass 
function. 2) Confinement and compositeness entail that 
the dressed-quarks comprising the bound state are not 
on shell and hence the two CP and T violating opera- 
tor structures that are indistinguishable in the free-quark 
limit yield materially different contributions to the EDM 
of the bound state. As a consequence these operator 
structures must be analysed and their strengths deter- 
mined independently in any model that provides for CP 
and T violation. Both these effects should be accounted 
for in using d n as a means of constraining extensions 
of the standard CKM model. These conclusions bear 



|0.85 h+ + 1.15 fq 



< 14.7 x 10" 



(42) 
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equally on the effects of operator structures we have ne- 
glected. 



APPENDIX A: IMPULSE APPROXIMATION 

As made clear in Sec. 2.3 of Ref. ]24] , when using an an- 
tisymmetrised product Ansatz for the nucleon's Fadde'ev 
amplitude the impulse approximation is 

5 

A,{q,P) = Al(q,P) + 2Y / Al(q,P), (Al) 

where: 

A^feP) - 3 J^L^e-l q; P)A 0+ (K) 

x^(£;P)A S ^p 3 + q,p 3 ), (A2) 

with K = £ + lp, p 3 = |p - £, Al{h,l 2 ) = 
x tr DF [A^(pi + g,pi)S(p2)] S(p 3 ) , (A3) 



the neutron EDM comes from A?, 
sum the contribution from the d-quark 



Of these five terms, the only u-quark contribution to 

the other four terms 
(The situation 
is reversed for the proton.) 

Our numerical results are obtained by evaluating these 
integrals using Monte-Carlo methods and the input spec- 
ified in Eqs. (0), ©, ©, ©, (J§) and ©. 
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where p x = -K + k, p 2 = -K - k, 6 = e clC2C3 e c 



is 



the pufcnrr contraction, and 



tt(pi,p 2 ,p 3 ) = A 0+ (p {12 })r +(ip[ 12] ;p {12} ) 

x^[ P{12} -2p 3 };P). (A4) 

A^ describes the photon probing the structure of the 
scalar diquark correlation, and contributes equally to 
both the proton and neutron. That contribution is triv- 
ially zero for the EDM. This merely reflects the fact that 
such a moment is forbidden to a scalar particle, so the 
only nontrivial contribution is that of the diquark's CP- 
and T-preserving electromagnetic form factor. The re- 
maining terms are 



X il(pi + q,p 3 ,p 2 ) S{p 2 ) {ir 2 y 
x A£'(pi,pi + g) (iT 2 )S( P3 ), 

x ft(pi,P2,P3) A^(p 2 + q,P2) 
x S( Pl ) S( P3 ) , 
d 4 k d 4 l 



Al(q,P) = 6 



— Q(pi,p 3 + q,p 2 ) 



(2tt) 4 (2tt) 4 

x tt(pi,p 2 ,p 3 ) S(p 2 ) S(pi) 
x A^(P3 + g,Pa)- 



(A5) 



(A6) 



(A7) 
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